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It is shown that Julia-Zee dyons do not admit slowly rotating excitations. This is achieved by 
investigating the complete set of stationary excitations which can give rise to non- vanishing angular 
momentum. The relevant zero modes are parametrized in a gauge invariant way and analyzed 
by means of a harmonic decomposition. Since general arguments show that the solutions to the 
linearized Bogomol'nyi equations cannot contribute to the angular momentum, the relevant modes 
are governed by a set of electric and a set of non self-dual magnetic perturbation equations. The 
absence of axial dipole deformations is also established. 



in 



O 

in 



X 



I. INTRODUCTION 



The main question addressed in this paper is whether Julia-Zee dyons admit rotational excitations. The investigation 
of this problem was motivated by some surprising results which we recently obtained for a class of self-gravitating 
non-Abelian soliton and black hole configurations. In jy we showed that the Bartnik-McKinnon solutions g] admit 
slowly rotating excitations. A two-parameter family of axisymmetric excitations of the static black hole solutions to 
the Einstein- Yang-Mills system was established as well. In addition to the charged, rotating black holes found in ||, 
there also exists a branch of uncharged, rotating black holes, as well as a branch of stationary - but not static - black 
holes with vanishing Komar angular momentum 

On the other hand, the situation was shown to be completely different in the presence of scalar fields Slowly 
rotating generalizations of (self-gravitating) solitons were excluded for a relatively large class of theories with non- 
Abelian gauge fields coupled to Higgs fields. In particular, the results obtained in [Q apply to the t'Hooft-Polyakov 
0^ ■ monopole and its self-gravitating generalizations. For black hole solutions of gauge theories with Higgs fields the 
""^ | situation is again different: Rotating excitations of static black holes generically exist; they are, however, necessarily 
Q" 1 charged. 

Since we are still lacking a deeper physical understanding of the facts mentioned above, we have been looking for 
(3Tj, other (not gravitating) examples which might help to find a clue. On the basis of our previous experience, we originally 
expected static dyon solutions to admit rotational excitations. The simplest examples are the Julia-Zee dyons, which 
are related to the Bogomol'nyi-Prasad-Sommerfield (BPS) monopole by a one-parameter family of hyperbolic rotations 
in internal space. 

The problem of small fluctuations around BPS monopoles has been examined some time ago by Mottola ||, Adler 
||, Weinberg ra], and completed in a comprehensive analysis by Akhoury et.al. The main emphasis was placed 
on the study of normalizable zero modes in the self-dual sector, because these are relevant to the structure of multi- 
monopole solutions. The moduli space of SU(2) monopole solutions carrying n units of magnetic charge was shown to 
be 4n-dimensional 0, ||. (See also pfj[ , |Tl| ], and (lj] for a generalization to arbitrary gauge groups and for further 
references.) As these studies are dealing with the self-dual sector, an investigation of the remaining zero modes seems 
to be necessary. This is also motivated by the following observations, which are obtained from general considerations. 

• The solutions to the linearized Bogomol'nyi equations - independently of whether or not they are physically 
acceptable - cannot give rise to a non-vanishing angular momentum. This is true for both BPS monopoles and 
Julia-Zee dyons. 

• The only excitations of BPS monopoles which can contribute to the angular momentum arise from perturbations 
of the time component, <5$ = 8A t , of the gauge potential; these will be called electric modes. 

• The only perturbations of Julia-Zee dyons which can contribute to the angular momentum are the electric ones 
and the non self-dual magnetic ones. (A mode will be called magnetic, if <5$ vanishes, and non self-dual, if it is 
a solution of the linearized field equations, but not of the linearized Bogomol'nyi equations.) 

The full problem, including the non self-dual fluctuations, was studied by Baake [|l3| in connection with the stability 
analysis of the t'Hooft-Polyakov monopole. In his work, Baake mainly focused on the negative fluctuation modes, the 
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absence of which he was able to prove by applying the Jacoby criterion. Since we are not aware of any other work 
devoted to non self-dual zero modes, we carry out a systematic, gauge invariant perturbation analysis in order to 
study the rotational excitations of BPS monopoles and Julia-Zee dyons. The emphasis in the present paper is mainly 
placed on the methods. The main result is, unfortunately, negative: Neither BPS monopoles nor Julia-Zee dyons 
admit slowly rotating excitations. 

A further motivation for studying non self-dual rotational excitations is provided by a theorem due to Taubes jl4) , 
according to which not every finite energy solution to the field equations in the BPS limit has to satisfy the first 
order Bogomol'nyi equations. Hence, the existence of physically acceptable excitations orthogonal to the Bogomol'nyi 
sector is not a priori excluded. However, the results of the present work imply that all non self-dual axisymmetric 
finite energy solutions, if they exist, are necessarily disconnected from the Julia-Zee dyons. This is, in fact, a week 
version of the original conjecture p5[ (the general form of which was disproved in Taubes' work |l4|]). It is, however, 
likely that configurations with unit winding number and discrete angular momenta exist. This is, for instance, the 
case for boson stars [ |l6j ]. 

This paper is organized as follows: In Sect. || we briefly review the symmetry which connects the PBS monopole 
solution with the one-parameter family of Julia-Zee dyons. In Sect. [II we show how to use this symmetry to reduce 
the perturbation analysis for Julia-Zee dyons to that for the PBS monopole. The main advantage of this consists in 
the fact that, after a hyperbolic rotation, the electric background field vanishes. This implies that - in the rotated 
system - the electric perturbations, <5<I>, do not couple to the magnetic ones. We then show that the non self-dual 
magnetic perturbations are governed by a system of first order equations for a one- form, SB. The latter comprises the 
perturbations of the Higgs field, SH, and the perturbation of the three-dimensional gauge potential, SA, in a gauge 
invariant way. 



In Sect. IV we present the decomposition of the gauge invariant perturbations Stf> and SB in terms of isospin 



harmonics. We also show that the expression for the angular momentum can be integrated, implying that only the 
boundary values of the perturbation amplitudes are relevant. The complete set of perturbation equations is derived 
in Sect. |y|. This consists of an even and an odd parity sector. Each sector comprises the electric equations for J<E>, 
the magnetic equations for SB (governing the non self-dual modes), and the inhomogeneous linearized Bogomol'nyi 
equations for SH and SA in terms of the source SB. 

In Sect. VI we discuss the odd parity perturbations and present the solutions of the complete set of equations in 
closed form. As the odd parity modes cannot contribute to the angular momentum, we conclude from the solutions 
that there exist no physically acceptable axial dipole deformations of Julia-Zee dyons. The more interesting even 



parity modes are discussed in Sect. VII. We show how to use the explicitly known solutions to reduce the magnetic 



problem to a standard Schrodinger equation. We also prove that the electric perturbations are governed by exactly 
the same equation. Since the latter has a non-negative potential, we are able to present a rigorous discussion of all 
modes. It turns out that there exist solutions (both electric and magnetic) which give rise to finite angular momentum. 
However, none of these modes are regular. 



II. BPS MONOPOLES AND JULIA-ZEE DYONS 

We consider stationary solutions to the SU(2) Yang-Mills-Higgs (YMH) equations with gauge potential and 
Higgs triplet H in the BPS limit (i.e., without Higgs self-interaction). The dimensionally reduced YMH action becomes 

S =\J {( F , p ) + (DH,DH) - (D$,D$) - [<S>, H} 2 } d 3 x, (1) 

where $ and A parametrize the electric and the magnetic components of the gauge potential, 

A (4) = $ dt + A. (2) 

The quantities F and D denote the field strength two-form and the gauge covariant derivative with respect to the 
three-dimensional magnetic potential A: 

F = dA + A A A, D$ = d$ + L4,$], DH = dH + [A, H). (3) 

(For arbitrary Lie algebra valued p-forms a the inner product is defined by (a, a)d 3 x = Tr{a A *a}, where * is the 
three-dimensional Hodge dual.) 

The perturbation analysis for Julia-Zee (JZ) dyons will be simplified considerably by the fact that the dimensionally 
reduced action (|l]) is invariant under hyperbolic rotations in the (H, $)-plane; that is, the transformation 
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H\ f cosh(7) sinh( 7 ) \ / H 
$ J ~~ * I sinh(7) cosh(7) J I $ 



(4) 



is a symmetry of the action (|l|) . 

In particular, a BPS monopole solution H = H mon , <E> = with magnetic charge P mon gives rise to a one-parameter 
family of JZ dyons, P = cosh(7)i/ mon , $ = sinh( 7 )i? mon , with magnetic charge P — cosh( 7 )P mon and electric charge 
Q = cosh(7) sinh(7)P mon . This is also seen from the field equations, 

*D*F= [$,D$] - [H,T)H], (5) 
*D*Dif =[$,[$, H]], (6) 

*D*D$ = [#,[$,#]], (7) 

which reduce to the monopole equations, D * F = — * [H mon , DiJ mon ] and D * DH mon = 0, for H — cosh(7)P mon and 
$ = sinh(7)iJ mon . 

It is worth recalling that the total energy is not invariant under the transformation (^). However, for fixed charges 
P and Q, defined by the flux integrals 

P = Jtt{HF} 7 Q = J Tr{H * D$} , (8) 

over the two-sphere at infinity, the energy assumes its global minimum for the corresponding JZ dyon solution. This is 
seen as follows: Using the field equations to express P and Q as volume integrals of Tr {DH A F} and Tr {DH A *D$}, 
respectively, the total energy may be expressed as follows [[17] , : 



E = - J {(F) 2 + (DP) 2 + (D$) 2 + [H, $] 2 } d : 



3 * 



= -[ \ (D$ - tanh( 7 )DP) 2 + (*F \ Bh] + [H, $] 2 1 d 3 a; 

2 J { V cosh (7) / J 

+ (Q sinh (^) + p ) . ( 9 ) 

cosh(7j 

where 7 is arbitrary and (F) 2 is a short-hand for (F, F), etc. From this one finds the bound (assuming, without loss 
of generality, that Q and P are non-negative) 



E > VQ 2 + P 2 = COsh 2 ( 7 )Pmon, (10) 

where equality holds if and only if A, H and $ are subject to the first order equation D<I>/ sinh(7) = DH/ cosh(7) = *F, 
which is exactly the Bogomol'nyi equation, 

* F = DH mon , (11) 
written in terms of the rotated fields H = cosh(7)_ff mon and $ = sinh(7)iJ mon . 

III. LINEAR PERTURBATIONS OF DYONS 

The perturbation analysis for the BPS monopole is simplified by the circumstance that the electric perturbation 
f^mon does not couple to the magnetic perturbations SH mon and i5j4 mon . This is an immediate consequence of the 
fact that the BPS background configuration is non-electric, $ m0 n = 0. 

Since the electric background field does not vanish for JZ dyons, the electric and the magnetic perturbations are 
coupled in this case. However, the linearity of the symmetry (|J) implies that all linear perturbations of JZ dyons 
can be obtained from the linear perturbations of the BPS monopole after a hyperbolic rotation with parameter 
sinh(7) = Q/P. It is, therefore, sufficient to consider the perturbation analysis of the BPS monopole. Before doing 
so, we compute the various contributions to the angular momentum. 
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A. Angular momentum 



The total angular momentum (along the symmetry axis) of a stationary YMH configuration is 

J= I T tip d 3 x, (12) 



where the relevant component of the stress-energy tensor in terms of the three-dimensional quantities is given by 

1. 



-Tr {[<!>, H]DH - *(D$A *F)} 



By virtue of the field equation (|) and the relations Tr {$[$, D$]} = and Tr {$[#, DH}} = Tr {[$, H]DH]}, we also 
find (after integrating by parts) 

T ttp = -±(*dTr{$*F}) v . (13) 

This shows that both the electric and the magnetic perturbations of JZ dyons contribute to the angular momentum, 
since 

ST tlp = -i(*dTr {,5$ * F + $ * SF}) V . 

(Note that the second term is absent if the electric background field vanishes, implying that only electric perturbations 
give rise to the angular momentum of a BPS monopole.) Since the dyon perturbations can be obtained from the 
monopole perturbations, we express the angular momentum in terms of the latter, using $ = sinh(7)_ff mon and 
5$ = smh(j)SH moa + cosh(7)6$ mon . With 

5T tv = cosh( 7 ) 5T? V + sinh( 7 ) 5T™\ (14) 

one finds 

$T t % = -\ (*dTr {6$ mon * F}) v , (15) 

ST^ = ~(5*dTt{H aum *F}) ip . (16) 

It is worthwhile noticing that both contributions to 5Tt v are separately gauge invariant. This is obvious for the 
electric part, since this is proportional to the perturbation of a field which vanishes on the background, namely $ mon - 
The same is true for the magnetic part, since the quantity dTr{i? mon * F} vanishes as well for a PBS background 
configuration. (Use *F = T)H mon to see this.) In fact, defining the one-form B according to 

B = DH moa - *F, (17) 

the magnetic contribution (|l^) to the angular momentum can be cast into the simple form 

8T^ = - l -{*^{5BA*F}) v , (18) 

which is manifestly gauge invariant, since, by definition, B vanishes for the BPS background configuration. 

The above expressions imply the following facts: First, the perturbation analysis for JZ dyons reduces to the 
perturbation analysis for BPS monopoles. Second, the electric and the magnetic perturbations of a BPS background 
contribute to the dyon angular momentum. Third, only the non self-dual modes, that is, the magnetic perturbations 
with SB 7^ contribute to the dyon angular momentum. 

The last statement reveals a fundamental difference between the perturbation theory of BPS monopoles and JZ 
dyons: Although the perturbation equations for JZ dyons can be reduced to the ones for the BPS monopole, the 
physical contents is quite different: While only electric perturbations can give rise to the angular momentum of a 
monopole configuration, magnetic perturbations need to be taken into account as well in the dyon case. Moreover, 
it is not sufficient to consider perturbations within the Bogomol'nyi sector, since the latter cannot contribute to the 
angular momentum of a dyon. 
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B. Linear perturbations of the BPS monopole 



Since the perturbation analysis of the JZ dyons can be reduced to the one for the BPS monopole, we shall now 
focus on the latter. In the following we omit the subscript "mon" indicating the monopole fields, that is, we write 
SH for SH mon , etc. Suppose that there is (at least) a one-parameter family of continuous deformations of the BPS 
monopole background, *F = DH, $ = 0. Then the tangent to this satisfies the linearized field equations. In order to 
linearize Eqs. (j^) and (jfy, it is very convenient to introduce the one- form field B defined in Eq. (H|). One may then 
write the first field equation in the form DB = D 2 H — D * F = [H, *DiJ — ^L - *D$], whereas the second field 
equation becomes D * B = D * DH — DF = *[$, [$, H]]. Hence, Eqs. (||) and (||) assume the form 

DB - [H, *B] = -[$, *D$], (19) 

and 

D* B = *[$,[$, H]\, (20) 

respectively. The linearization of the field equations (^), (pJ|), fl20| ) is completely trivial, since both the electric field «3> 
and the magnetic one-form B = DH — *F vanish for a BPS background. Hence, the linearized field equations involve 
only the gauge invariant perturbations (5$ and SB: One immediately finds the result 

electric perturbations: D * DS& — *[H, H]], (21) 
magnetic perturbations: D8B = [H, *5B], D * SB = 0, (22) 



where SB is obtained from the definition (|17|), that is 

SB = DSH - *DSA - [H, 5 A}. (23) 

(Here and in the following all quantities without a "5" refer to background fields.) Before we consider the harmonic 
analysis of Eqs. (pl|)-(p3|), we note the following: 



The linearization of the BogomoPnyi equation (|ll|), SB = 0, has been studied extensively in the literature. 
The solutions to SB = are, however, only a subset of the general magnetic perturbations. The full magnetic 
perturbations are governed by the second order equations for 5 A and SH, which are equivalent to the first oder 



equations (22) for SB and the inhomogeneous equation (|23|). In particular, we have already argued above that 
only the non-trivial solutions SB ^ to Eq. (|2^) can contribute to the angular momentum [see Eq. (p"8|)]. 

In order to find the general magnetic perturbations, one proceeds in two steps: First, one has to solve the system 
( f22| ) for SB. Once SB is known, it remains to solve the inhomogeneous linearized Bogomol'nyi equation (|2^) for 
8 A and SH. This is achieved by using Green's method, also taking advantage of the explicitly known solutions 
to the homogeneous equations, SB = 0, derived in || and ||. 

Since the background BPS configuration has vanishing B, the magnetic perturbation SB is manifestly gauge 
invariant. This is also verified by using the general behavior of the perturbations 5 A and SH under gauge 
transformations generated by a Lie algebra valued scalar field %: 

5 A -> 5 A + B% SH -> SH + [H, x] ■ (24) 

Hence, SDH -> SDH + [DH, X ], and SF -> SF + [F, X ], implying that SB -> SB + [B, X ] = SB. 

The second equation in (E2p is a consistency condition for the first one: Indeed, applying D on the first equation 
and using D 2 SB = [F, SB] on the LHS, and [DH, *SB] = [*DH,SB] = [F,SB] on the RHS, yields the necessary 
condition [H, D * SB] = 0. 



IV. HARMONIC ANALYSIS 



Since the unperturbed BPS solution is spherically symmetric, we perform a multipole decomposition and rewrite 
the electric perturbation equations ( pl| ) and the magnetic ones (^2|), (|23|) as systems of ordinary differential equations 
with respect to the radial coordinate. Using these equations, we show that the angular momentum integral can be 
computed exactly. Hence, the total angular momentum arising from electric and magnetic perturbations is determined 
by the asymptotic behavior of the gauge invariant amplitudes S$ and SB, respectively. 
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A. Isospin harmonics 



The basic fields H, <I>, A, the auxiliary held B, and their perturbations, are functions and one-forms with values in 
the Lie algebra su(2) of the gauge group SU(2). Let us start by considering such functions on the two-sphere S 2 . A 
convenient basis, reducing the natural representation of SU(2), is obtained by taking the inner product of the vector 
spherical harmonics Yj M with the basis r = <r/(2i) of su(2) (where er are the Pauli matrices) 



JM 



The isospin harmonics Cj M have total angular momentum J and fixed parity (— 1) L . Instead of the Y 



usual to consider the basis Y^ M (with A = 0, ±1). For A = and A = 1 these vector harmonics are transverse, while 
they are longitudinal for A = — 1 (with respect to the radial unit direction f) 



JM 



(25) 
it is also 



Yj° M are also called electric and magnetic multipoles, respectively. They are obtained by applying certain differential 
operators on the ordinary spherical harmonics Yjm, while the longitudinal harmonics are given by Y 



The transverse harmonics Yj M and 



(-i) 

JM 



= fY 



JM 



(see, e.g., or ptj). The formulas for the Yj M can readily be translated into the corresponding formulas for the 



isospin harmonics Cj X M = r • Y^ M (with A = 0, ±1). One finds 



ni- 1 ) — v 



c 



JM 
(0) 



c 



JM 



( + 1) 
JM 



i 



= (dr r , JdF 



JM 



(26) 



(dr r , dY JM ), 



where r r = r ■ r. Here ( , ) and * denote the inner product and the Hodge dual with respect to the standard metric 
on S 2 . (Also note that the spherical components of the r obey the relations dr r = r^d'd + t v sim?d</? and [t§, r„] = r r ; 



see Appendix ^.) In terms of the isospin harmonics, the well-known relations between the vector harmonics Y 



and Yy£j become 



L 

JM 



^JM 



r j 

^JM 



1 



V2J+1 



C 



C 



7-1 
JM 



(0) 
JM' 

1 



V2J+ 1 



(-1) 

JM 



(27) 



By construction, the isospin harmonics Cj M , Cj^ are eigenfunctions of the spherical Laplacian, A = Sdid, and of 
the parity operator, P: 



kC L JM = -L{L + l)C L 



JM, 



(28) 



PC 



JM 



(-1) 



1) L C L 



JM, 



(29) 



(The exterior derivatives of the isospin harmonics Cj M and Cj^ and their S 2 duals are 



where L = J, J ± 1 

particularly convenient for analyzing perturbations of Lie algebra valued one-forms 
of monopole harmonics we refer to p2|.) 



For the general theory 



B. Perturbation Amplitudes 

Since rotational modes are our primary concern in this paper, we now focus on the sector J = 1. For the c[q 
(A = 0,±1) we use (with some change of normalization) the letters X, Y and Z. A convenient basis of J = 1 isospin 
harmonics then is 

X = tv K, where K = cos$, 
V2 Y = (dr r , dK) = -T,? sim?, 

V2 Z = - (dT r , *dK) = t v sim?, (30) 



G 



where X and Y span the even parity sector, while Z has odd parity. The su(2) valued electric perturbation function 
can, therefore, be expanded as = (5$ cvcn + 5$ odd , with 

(5$ cvcn = - {ct>-X + 0+Y) , 
r 

(The factor 1/r is introduced for convenience; see, e.g. Eqs. (f45|), (Efj). Throughout this paper, all amplitudes 
furnished with a tilde refer to the odd parity sector, which is relevant for deformations.) A similar expansion holds 
for SH; however, unlike <5$, SH is not gauge invariant; see Sect. VC and Appendix ^. 

Turning to Lie algebra valued one-forms, we note that the exterior derivatives of the basis functions X, Y and Z 
can be expressed in terms of the derivatives of r r and K = cosi? oc Y\$. (This is a peculiarity of the J = 1 harmonics, 

for which dCflj = (V2dY + dX)/VS = 0.) One finds 

dX = - V2 dY = 7> dK + K dr r , 

*V2dZ = T r dK-KdT r . (32) 

As the parity operation commutes with the exterior differentiations and anti-commutes with the Hodge dual, one can 
expand the su(2) valued magnetic perturbation one-form SB as SB = SB evcn + SB odd , with 

5B even = 4r (b-X + b+Y) dr + fa T r dK + f3 2 Kdr r , 

SB odd = \(bZ)dr + fa ir r dK + fa *Kdr r , (33) 

where b, b±, fa^. and fa$ depend on the radial coordinate r. (Again, a similar formula holds for 5 A. In contrast to 
SB, 8 A is not gauge i nvar iant, implying that not all coefficients in the expansion of 5 A correspond to physical degrees 
of freedom; see Sect. VC and Appendix [5]) . 

At this point we also recall that the background gauge potential and Higgs field are parametrized in terms of two 
radial functions w(r) and h(r) (see Appendix^), 

A = [1 - w(r)] *dr r , H = h{r)r r . (34) 

Since r r is an eigenfunction of the spherical Laplacian, dSdi> = — 2T r dQ, the background field strength becomes 
F = -dw A id-n r + (w 2 - l)r r dfi. The BPS equations, F = *DH, thus read 

w' = wh, r 2 h' = w 2 -l, (35) 

with the globally regular solution 

, . r , . . 1 cosh(j') . . 

wr ) = ^iTT' h ( r = r rH- 36 

smn(r) r sinii(r) 

For later use we also note that the second order equation for h can be integrated, which yields the useful relation 

h'=h 2 - — -l. (37) 
r 



C. Integration of angular momentum 

We now show that the total angular momentum SJ = cosh(-f)5 J cl + sinh(7)<5J mg can be expressed in terms of the 
values of the gauge invariant perturbations <5$ and SB at the origin and at infinity. According to Eqs. (15) and (|Tg|), 
the electric and the magnetic perturbations give rise to 

SJ cl = ~\! (*dTr {(5$ ovcn * F}) v d 3 x (38) 
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and 



SJ ms 



(*Tr{<5B ovon A *F}U 3 i, 



(39) 



respectively. Here we have already used the fact that only the even parity sector contributes to the total angular 
momentum. In order to express the above integrands in terms of the radial amplitudes cf>±, b± and /3i,2, we first note 
that the background field strength can be written in the simple form 



* F = w'dr r + h'r r dr . 



(40) 



Taking advantage of the trace formulas Tr {Xr r } = -K/2, Tr {Yr r } = Tr {Zr r } = 0, and Tr {Xdr r } = 0, Tr {Ydr r } 
*Tr {Zdr r } = — dK/{2 y/2), it is now not difficult to compute the above integrands from the expansions ( |3l|) and (p^ 
One finds 



*dTr{5$ cvcn *F} 



1 




r 2 ( w'4 




UK, 




V2\ r 




1 

2T 2 " 


r 2 h'(3! - 


1 " 

T2 wb \ 


*dK. 



(41) 
(42) 



With K = cost? we have *dK = — sin 2 i?d<yj, which shows that the above formulas yield the <p components appearing 
in the integrands of Eqs. ( p8| ) and (|3^). It is an interesting fact that the above brackets can be written as radial 
derivatives. This enables one to perform the angular momentum integrals and to express 5J el and <5J mg in terms of 
the values of ^$ cvcn and SB cvcn at the origin and at infinity. In order to see this, one has to use the perturbation 
equations in the harmonic decomposition, as given in the next section. Considering the magnetic part, one uses the 

which enables one to eliminate (3i in Eq. (|42"1). Also 

.21 



first two equations in (^9| ) to obtain 2(1 — w 2 )j3\ = b'_ + \/2wb'^ 



taking advantage of the background equation (|35|), one then has \r 2 h! B\ — w'b+/y/2] = — [b- + y/2wb + ]' /2. A similar, 
but more complicated manipulation uses the second order equations (|4^) to write the electric contribution (|4l| ) in the 
desired form; see Appendix [F]. The two contributions (|3^) and ( |39| ) to the angular momentum finally become 



5J cl = - 



6J n 



(l-w 2 - 2rh) + r 2 h(f)'_ + V2wrh(f)+ 
6_ + V2wh 



(43) 
(44) 



V. PERTURBATION EQUATIONS 

Using the expansions (|3l| ) and (|33|), as well as the tools developed in Appendix ^| and|c[ it is now straightforward 
to write down the system of perturbation equations. This consists of Eq. (|2l|) f or the electric perturbations <5<I>, Eqs. 
(^||) for the magnetic perturbations SB, and the inhomogeneous BPS Eqs (p3|) for SH and SA. 



A. Electric perturbations 

For the electric perturbations ([H]), governed by Eq. (|2l]), one finds the differential equations 

(45) 



1_ ( 2(w 2 + 1) -2V2w 
r 2 V -2^/2w (w 2 + 1 + r 2 h 2 ) ) \ </>+ 



in the even parity sector, and 



in the odd parity sector. Here we have used Eq. ( |Cl| ) to compute the LHS of Eq. (|2l|), and [r r , X] — 0, [r r , Y] = —Z, 
[T r , Z] = Y to obtain the RHS: [H, [5$, H]] = tfr^^+Y + %Z). 



I" = ^(w 2 + l + r 2 h 2 )4> (46) 
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B. Magnetic perturbations: SB equations 



In order to determine the magnetic perturbations SB, we first write the decomposition ( |33| ) in the form 

5B = -^bdr + B, (47) 

where the one-form B is tangential to S 2 . In terms of b and B, the magnetic perturbation equations ( p2|) assume the 
form 

[H, b] = *£)£, 
B' = r- 2 t>b - [H,*B], 

b' = -*I)iB. (48) 

Here we have used the fact that the unperturbed gauge potential has no radial component, implying the decomposition 
D = dr A d 7 ■ + D for the covariant derivative (see Appendix |b|) . Taking advantage of the formulas given in Appendix 
it is now not hard to obtain the sets of differential equations for the radial functions parametrizing 5B even and 
SB odd . One finds 

b'_ = 2(01 + wfc), b' + = -V2(wpi + ft), 



II) I 

r 2 (3' 1 = b---^=b + , r 2 (f3 2 + hf] 2 ) = wb--^=b + , 



for the even parity sector, and 



hb + = V2(p 2 - wfr) (49) 

V = V2{w0t - /3 2 ), 
r 2 /3i = ^U, r 2 0' 2 + hp 2 ) = ~b, 

= /?i+u;/32, hb= V2{wh +/3 2 ) (50) 

for the odd parity sector. We note that both sectors contain constraint equations, reflecting the fact that the second 
equation in ( |2^ ) is an integrability condition for the first one. 

At this point we also note the following, somewhat surprising fact: The scalar magnetic amplitudes b± and b satisfy 



the same set of second order equations (|45j), (46) as the electric amplitudes cj>± and (f>, 

-2V2w (w 2 + l + r 2 h 2 ) J \b+ J ' 



b" = —(w 2 + l + r 2 h 2 )b. 



(52) 



This follows from the arguments given in Appendix H and is also verified directly from the above equations. Using 
the odd parity equations J5(i| ) we have, for instance b" = V2(w'/3i + w/3[ — f3' 2 ) = [h 2 + r~ 2 (w 2 + l)]b, where we have 
used the Bogomol'nyi equations ( p5| ) for the background fields h and w. 

We also point out that not all solutions to the second order equations ( |5l| ) and ([32]) satisfy the first order equations 
( p9| ) and (50). In fact, it is not hard to see that the solution spaces defined by Eqs. ( f49| ) and (50) are three- and 
onc-dimensional, respectively, rather than four- and two-dimensional. 
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C. Magnetic perturbations: Inhomogeneous BPS equations 



In order to write out the inhomogeneous Bogomol'nyi equations (23), we need the harmonic decomposition of the 
fields SH and SA. The fact that the latter are not gauge invariant enables us to get rid of certain amplitudes. In 
Appendix ^ it is shown that - up to a pure gauge - the harmonic decompositions of SH and 5 A assume the form 



SH even 
SH odd 



1-X 



■J+Y, 



(53) 



and 



SA even = ai *T r dK + a 2 iKdT r 
SA odd = &i T r dK + a 2 Kdr r , 



(54) 

respectively. The radial functions j±, 7, a\ y 2 and aj.,2 are gauge invariant, up to a one-dimensional set of residual 
gauge transformations in the even parity sector, 



7- 
a.1 



Ct2 



7+ 
w 

1 

V2 



7+ + hc 3 , 



C3, 



C3, 



(55) 



and a two-dimensional set of residual gauge transformations in the odd parity sector, 

7 — > 7 - hC2, 



a 1 



OL2 



Oi\ + Cl 



V2 



(-2: 



012 



WCl -= C 2 , 

v2 



(56) 



where c\ , C2 and C3 are arbitrary constants parametrizing the residual gauge freedom (see Appendix [d]) . In terms of 
the gauge invariant source terms b± 7 b, (3\^ 2 and /?i.2, and the (almost) gauge invariant amplitudes introduced above, 
the inhomogeneous linearized Bogomol'nyi equations (|23| ) eventually become 

r 2 j'_ + 2{ai + wa 2 ) 



r 2 i + /V2 - {wai + a 2 ) = b + /V2, 
a[ + 7- - wj + /V2 = 0x, 
W7- -7+/V^= P2, 



a 2 — /ia2 



(57) 



in the even parity sector, and 



r 7 



(a 2 - wi. 



6/V2, 



— 6e 1 + wj/v2 = Pi, 
\' 2 + ha 2 — 7/V2 = fo, 



(58) 



in the sector with odd parity. 

For vanishing RHS, the above equations are the linearized Bogomol'nyi equations, which have been studied in the 
literature. Using the background equations (|3^), it is easy to verify that the residual gauge mode 7_ = 0, 7+ = y/2h, 
a\ = w, oi2 = — 1 satisfies the homogeneous equations ([57]), while the residual gauge modes 7 = y/2h, a\ = w, &2 = 1 
and 7 = 0, <5i = 1, a.2 — w are solutions to the homogeneous equations (58). 



VI. ODD PARITY MODES 



We shall now solve the perturbation equations. We start with the odd parity sector, for which all solutions can be 
obtained in closed form. We emphasize, however, that t his se ctor is of minor importance, since the odd parity modes 
cannot contribute to the angular momentum. In Sect. VIA we compute the magnetic amplitudes SB odd , which we 
use in Sect. VI B as source terms to obtain the perturbations SH odd and i5A odd . In Sect. VIC we finally compute the 



electric perturbations (5$° 



da 
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A. Solutions to the SB equations 



In order to compute the source term i5i? odd , we have to solve Eqs. (50) for the amplitudes b and j3\^ defined 
in Eq. (|3^). Using the last two equations in ( |50| ) to express f3i and 02 in terms of b, the first equation becomes 
b'/b = h(w + w^ 1 )/(w — w^ 1 ), which is trivial to solve, since the numerator is the derivative of the denominator. 
Hence, the only solution to Eqs. ( |5C| ) is 



-, V2Pi=h, V2j3-2 = — . (59) 



h 

w 

w w 

Inserting this back into the expansion (|33|), and using the background equation (^5|) for h! and the formula flCl| ) for 
DZ, yields the simple result 

SB odd = -D(hZ). (60) 
w 



B. Solutions to the inhomogeneous BPS equations 



Now that the source terms for the linearized Bogomol'nyi equations ( p8| ) are known, we can proceed and solve 
the inhomogeneous problem. Since the homogeneous equations admit three solutions, two of which are the residual 
gauge modes 7 = V2h, &% = w, 0.2 = 1 and 7 = Q, u\ = 1, 0.2 = w, we need to find the remaining solution of the 
homogeneous problem and a solution of the inhomogeneous equations. This is achieved by deriving a third order 
equation for 7. In fact, since 7 = is a residual gauge mode of Eqs. (p8|), the differential equation for 7 will be of 
second, rather than third order. Moreover, using the second residual gauge mode, one eventually ends up with a first 
order equation. First, one easily finds from Eqs. (K8h 



(r 2 f) - h (r 2 f) - (w 2 + 1)7 



-hb. 



(61) 



where \/2(w(3i— fc) = b' was used on the RHS. Now using the second residual gauge mode 7 gauge 
part of the above equation can be cast into the following first order equation for (7//1)': 



h, the homogeneous 



2„2 



K 2 r 



0, 



with the solution 7 cx h J w/(rh) 2 . The integration can be performed by using the relation [w/(r 2 h)]' — w/{rh) 2 , 
following from the background equation (37). Hence, the only non-gauge mode of the homogeneous equations (pq) is 



~hom _ ^/2a\ om = h--, V2<5^° 



w 
r 



(62) 



(In order to verify that this solves the homogeneous part of Eqs. (|58|), one uses again the first order equation Q3^ 
for the background field h.) We may finally use the two solutions 7^°™ and ^saugc _ ^ ^ so i ve th e inhomogeneous 
equation ( [31] ) with source term X — —hb — —h(w — w^ 1 ): 



~inh 



*r ~ 2 (^(i)7 sausc + ^( 2 )7 h ° m ) 



(63) 



with = j hom /W and (ip) = -7S au g°/IU, where W = ^uge^homy _ ^hom^gaugey is the Wronskian of the two 
homogeneous solutions. A short computation yields W — w/r 2 , and hence fin) = 1, (i( 2 ) = —r 2 h/w. We thus end up 
with 



~inh = ^_ 



r 2 h hti 



hti 



(64) 



This shows that the physical modes describing magnetic perturbations with J = 1 and odd parity form a two 
parameter family. In particular, the perturbations of the Higgs field become 
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(The arbitrary constant C2 reflects the fact that the source terms SB are themselves solutions to a homogeneous 
set of equations, implying that the inhomogeneity in Eq. ( |6l| ) is only fixed up to a multiplicative constant.) Since 
the self-dual solution ^ hom diverges like 1/r 2 near the origin, while the non self-dual part 7 inh diverges like J e r jr 
at infinity, we conclude that there exist no small magnetic perturbations of BPS monopoles and JZ dyons with odd 
parity. 



C. Solutions to the 5$ equations 

The electric perturbations 4> with odd parity are governed by Eq. (^). Since the magnetic amplitude b fulfills the 
same second order equation, we immediately conclude from the solution (|5^) that 

= w (66) 

w 

solves Eq. (^). [In fact, using (w ± w^ 1 )' = h(w =F w^ 1 ), one has (w — w^ 1 )" — h'(w + w^ 1 ) + h 2 (w — w^ 1 ) — 
[{w 2 + l)/r 2 + h 2 ](w — W 1 ).] The second solution is given by (j)^ = 4>^ J[4>^]~ 2 dr. The integral can be carried 
out, and yields 

4,™ =±L+r\--- (67) 

r \ w J w 

[Using the background equations ([35j) it is not hard to verify that this is indeed the second solution to Eq. (E3).] The 
electric perturbations with odd parity are, therefore, 

5$ odd = - (Ci0 (1) + C 2 (2) ) , (68) 

which remains finite for r — > 00 only if C\ = C2. However, as (tp 1 -* + <p 2 ^)/r diverges like 1/r 2 in the vicinity of the 
origin, we conclude that there exist no small electric perturbations of BPS monopoles and JZ dyons with odd parity. 



VII. EVEN PARITY MODES 

A. Solutions to the SB equations 

In order to solve Eqs. (|4^), we first note that the equation for b' + is a consequence of the remaining ones. Eliminating 
b + by using the last equation in (^) , we obtain a system of three first order equations for 6_ , (3\ and . It is then 
straightforward to decouple these equations, which yields a third order equation for b~. Since b- enters this equation 
only via its derivatives, one concludes that 6_ = const, is a solution. In fact, one easily verifies that (any constant 
times) 

6 (0 > = 2, &f = V2 (w + IV p{ 0) = -h, ^ = A (69) 
\ w J w 

solves Eqs. (^). In order to find the remaining two solutions, it is convenient to define the quantities 

E = hb- + (fa - wfe) , 

A = h~ x (/3i + wfo) • (70) 

Since S and A vanish for the solution (^), it is possible to derive a system of two first order equations for these 
quantities. Using the three equations for (3\ and 02, one finds after some manipulations S' = 2h 2 A and A' = 
h 2 r 2 (w 2 + 1) _1 E, which also yields the following second order equation for E: 

h! w 2 4- 1 

S" - - 2^^£ = 0. (71) 

h r 

Considering the quantity E//i, we obtain a Schrodinger equation with potential P(r), 
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w 2 + 1 (h^ 1 )" 

(S/ft)"=P(r)(S/ft), whereP(r)=2 ±- + (72) 

r A h 1 

Having solved this equation, one obtains the magnetic amplitude &_ from the definition of E and the first equation 
in (49). In order to find b + one uses the last equation in ( [49| ) and solves Eqs. ( fTOj ) for /3i and /3 2 . This yields 



f E' 

^-H" + =) (M - E) -5K"-=) !! '- <73) 

The formulas for /3i and /?2 in terms of E arc n ot given here, since SB can be expressed in terms of 6_ and 6+ alone. 
This is seen as follows: Using the relations (|ci"|), the terms tangential to S 2 in the expansion (33) for SB e 
written in the form 

fa T r dK + (3 2 KdT r = 



w 2 — 1 

1 



(hb- -E)DX + (hb + )f)Y 



w 2 — 1 

where we have used the last equation in ([l9]) and the definition ( f70| ) to get rid of f3\ and fa. Now using Eq. ( f73| ) for 
&_, we have E' = and thus (/i&_ — E)' = h'b-, which enables us to write the term proportional to X in Eq. ( |3^ ) 
as r~ 2 b-dr = (w 2 — l)~ 1 h'b-dr = (w 2 — l)~ 1 d(/i&_ — E). Hence, the terms proportional to X and DX combine to 
an exact covariant derivative, which finally yields the result 

SB cvcn = — {D [{hb- - E) X] + b+B [hY]} . (74) 



This shows that all three magnetic modes with even parity are obtained from Eqs. (|72j), (|73|) and Q74|). In particular, 
the trivial solution E^ = of Eq. (172) gives rise to the solution (p9), which yields 



SB^M = | D {hX) + T2 ( W + -)D (hY)j . (75) 

The two non-trivial solutions of Eq. ( |72| ) are not (yet) known in closed form. However, their qualitative behavior 
can be discussed rigorously: The potential P(r) is positive definite for all finite values of r. (In fact, using the 
background equations to compute (h^ 1 )" , one finds from Eq. d73) 



1 fh'\ 2 (\ h'^ 2 



P W = ^ + UJ + {r + ^) ' (?6) 
which is manifestly non-negative, and vanishes only for r — * oo.) As h'/h = 1/r — 2r/15 + C(r 3 ) in the vicinity 



2 



of the origin, we have P(r) = 6/r 2 + 0(1), implying that the fundamental solutions behave like E^ /h oc 1/r 
and E( 2 V^i oc r 3 . For r oo one has P(r) = 2/r 2 + 0{l/r 3 ), which yields T,^ /h oc 1/r and Z^/h oc r 2 . The 
monotonicity property of E//i, following from the positivity of the potential ((7^), enables one to conclude that the 
solution which diverges at the origin remains finite at infinity, and vice-versa. (Also we have used global existence, 
following from the linearity of Eq. (|7^) and from the finiteness of the potential for r ^= 0.). Since h behaves like r 
near the origin and approaches the constant value —1 at infinity, the two non-trivial solutions of Eq. ((72|) behave as 



follows: 

E (1) ocr-\ E (2) ocr 4 , for r -»■ 0, (77) 

and 

E (1) ocr-\ E (2) ocr 2 , for r -> oo. (78) 



In Sect. IV C we have shown that the angular momentum can be expressed as a boundary integral. The relevant 
quantity appearing in Eq. (Q) is the difference of (6_ + \/2wb + ) between infinity and zero. For the solution E^ this 
quantity remains finite at infinity, whereas it diverges logarithmically at the origin. (Note that the leading power, 
1/r 2 , cancels in the above combination.) On the other hand, (fo_ + \/2wb+) remains bounded at the origin for the 
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solution £( 2 \ whereas it obviously diverges like r 2 at infinity. Hence, neither of the two non-trivial solutions to 



Eq. (|72|) gives rise to a finite angular momentum. The fact that we are able to decide this without solving the 
inhomogeneous equations (^) for SH and 8 A follows from the observation that the angular momentum depends on 
the magnetic perturbations only via the field SB. Moreover, only the boundary values of the gauge invariant quantity 
SB are needed to obtain the magnetic contribution to the total angular momentum. 

Surprisingly enough, the third solution, given in closed form in Eq. (|69|), does give rise to a finite angular momentum, 
although the amplitude b+ diverges at infinity. By virtue of Eq. we have 



gjmg = 

3 



b_ + V2wb _, 



2tt 



(79) 



In order to decide whether this is an acceptable perturbation, we have to compute the physical fields SH and 8 A. 



B. Solutions to the inhomogeneous BPS equations 



We now discuss the inhomogeneous linearized Bogomol'nyi equations (|57|). These can be written as a fourth order 
equation for either of the four variables 7±, 0:1,2, parametrizing SH cvcn and SA cvcn . Since 7_ vanishes for the residual 
gauge mode of the system ( |57|) , the equation for 7_ is only of third, rather than fourth order. One finds 

(80) 



(r 2 7 ^)" - 2h (r 2 7 l)' - 2 (l + 2w 2 ) i_ + Ah (l - w 2 ) 7 _ = 1[SB C 



where the inhomogeneity is an expression in terms of b± and /3i,2- Using Eqs. ( f49|) for these amplitudes, it is possible 
to express T[SB° vcn ] in terms of b + and alone: 



I[SB C 



2w 



w{hbJ)' + V2(hb+)' 



(81) 



The solutions to the homogeneous problem, that is, the solutions to the linearized Bogomol'nyi equations are known 
in closed form ||, ||, In fact, they can be expressed in terms of the quantities w and h. Using the background 
relations (35), it is not difficult - and not particularly pleasant either - to verify that the three solutions of the 
homogeneous Eq. ( |80| ) are 



„(l) -h' - 



(2) 



(3) 

7- 



= h< = -( w 2 -l), 

4K-(rk-l)] 
2 (rh - 1) . 



(82) 



Concerning these solutions of the homogeneous linearized PBS equations, we note the following. 



(2) . . ... (3) . . . (2) 

• The amplitude 7_ is not finite at the origin, while 7_ becomes unbounded at infinity. Hence, neither 7_ nor 
7_ give rise to small perturbations of SH cvcn . 

• The third solution, 7_ ' , does give rise to an acceptable physical mode. The latter corresponds to a translation 
along the z-axis. This is seen by differentiating the background field H — r r h with respect to d z = cos{)d r — 
r^ 1 sm$d#, which yields d z H = h'X + r^ 1 h^/2Y, and hence 7_ = hf. (Note that the coefficient in front of Y is 
not the amplitude 7 + introduced in Eq. (p3|), since the latter was defined in a gauge where SA is tangential to 
S 2 ; see Eq. (]54). The only quantity which is not affected by the corresponding gauge transformation is 7_; see 
also Appendix ^.) 

• The remaining two physical zero modes in the sector J = 1 correspond to translations along the x- and y-axis. 
They do not occur in the above calculation, since, for reasons of symmetry, we have restricted the harmonic 
decompositions to the magnetic quantum number M = 0. 

• We recall that none of the above solutions to the linearized BPS equations can contribute to the angular 
momentum, since only the field SB, describing the non self-dual perturbations, enters the expression (p~q). 
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It remains to find the solutions to the inhomogeneous equation (|0|). Since all solutions to the homogeneous problem 
are known, we can apply standard methods to obtain the particular solution j_ ■ For given inhomogeneity I one has 

7 1 l nh = E^ ) /MW Jdr, (83) 

(k) 

where the three quantities fi^k) are obtained from the homogeneous solutions 71 by 

^ ~ (m), (n) v , (<)y, • ^ M > 

EmnlT- (7- ) (7- ) 

A rather lengthy computation yields the value — 24w 2 /r 4 for the Wronskian in the denominator, and then 



A*(i) = ^ (l + r 2 + sinh 2 (Y)) , /U (2 ) = ^ ^sinh(r) cosh(r) - r - ^r 3 ^ 



M(3) = -^. (85) 
In the previous section we have argued that only the solution (|69|) of the equations (|49|) for SB cvcn gives rise to a finite 



angular momentum. Hence, it remains to compute 7^" h for the source term given by b_ — 2 and b + — y/2(w + w 1 ) 
Using the expression (|8l|), we immediately have 



J[5B cvca ^] = -4 ( h 2 + 2w2 + - ) . (86) 



The solution 7" lh is now obtained from Eqs. (p3|) , (p4|), ( p5|) and (p6|) . An expansion in powers of r reveals that Y- h 
diverges like 1/r in the vicinity of the origin. Since 7^ diverges like 1/r 3 , while 7^ and are well-behaved for 
r — > 0, there is no combination of 7™ 11 with a homogeneous solution ( p^ ) which remains bounded at the origin. Hence, 
we conclude that there exist no bounded magnetic perturbations 5H , 6 A, which give rise to finite angular momentum. 



C. Solutions to the <5<E> equations 

It remains to discuss the electric perturbations with even parity. The latter are governed by Eqs. (^) for 0_ and 
4>+. We have already argued that three solutions of these equations coincide with the magnetic solutions b±\ b±* 
and b±\ discussed in Sect. |VII A . In order to find the remaining solution, it is convenient to write the system ( |45| ) 



in the form of a second order equation with two inhomogeneities. The manipulations by which this can be achieved 
are discussed in Appendix [|]. Introducing the quantity S in the same way as in the magnetic case [see Eq. (|73|)], 



E' = h<f>'_, (87) 



one eventually finds the equation 



h' ~ . w 2 + 1 - . w 2 + 1 
S" - 2-S' - 2 — E = -ka ti - k 2 —, 

h r -2 r z 

where fco and ^3 are integration constants. The homogeneous part of this equation coincides with the corresponding 
magnetic equation (fn]). 

The particular solution for k$ — is E = koY>(°\ with E' ) = 1. This yields (f>^ = const., which coincides with 
the magnetic solution ([39|). The two solutions to the homogeneous problem, E^ and E^ 2 ', say, coincide with the two 
remaining magnetic solutions obtained from the homogeneous equation (|7l|). The additional solution, E — fc3S^ 3 \ 
which is not present in the magnetic case, is the particular solution for ko = 0. Again, this can be given in closed 
form by introducing the quantity S = E//1 2 . Using again the background equation d37j ) for ft,, a short computation 
shows that the LHS of Eq. ( |88|) assumes the form (S'h 2 )' + 2h'S. Hence, the desired solution is S = —k^/2, that is, 
£(3) = -h 2 /2. By virtue of Eq. @ this yields ^ oc h. 



We thus conclude that the four electric perturbations with even parity are given by 
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do) 



V2 



= 2h, 



V2 



(rh)'. 



,(1.2) _ 



V2 



w H 



(1:2) 



„2/i,(l,2) v , 



2w 



(89) 
(90) 
(91) 



where & 1<2 > are the two non-trivial solutions to E" — 2(/i'//i)£' — 2r 2 (w 2 



1)E = 0. The angular momentum is 



obtained from Eq. ([13|) and the above solutions. One finds 
6J**> = ^.[rh{h-rh')]Z, 



5 jel(l,2) = _: 



b! 



J h I 'J_,..,, 1^^.(1.2) 



-(rh- 1) 



£' 2) )'-^' 2) )" 



(92) 



It is easy to see that the only combination of <fr( ' and (j)^ which gives rise to finite angular momentum is their 
sum. However, since the amplitudes entering <5$ are 4>+/r and 4>~/r [see Eq. (|3l|)], the perturbation <5$ obtained 
from (f)^ + (f)^ diverges at the origin like 1/r. 

The behavior of f)W given in Eqs. @ and (|7|) implies that ^ = C(l/r 2 ) as r -> 0, and = 0(l/r) as 
r — > oo. Using this in the above expression shows that the angular momentum is again finite, | <5J el W |< oo. As 
above, the perturbation (5$ is, however, not bounded at the origin. 

The solution <)P^ diverges like r 2 as r — > oo. However, the leading order terms in the expression for <5J c1 ' 2 ' cancel, 
and so do the next-to-leading order terms. Hence, like <5J el (°) and SJ 01 ^, SJ cl ^ diverges only with the first power of 



r, implying that there exist linear combinations of (p^ with 4> { " 1 (or (fif) which give rise to finite angular momentum. 

Since <f>_ /r is not bounded at the origin, while (j>^ /r is bounded, only linear combinations of <f>^ with <fi^ need to 

be considered. However, the latter give rise to perturbations (5$ which are not bounded at infinity. (Note that <p\ /r 

behaves like e r /r, whereas 0+V r grows like e r .) 

We therefor conclude that all electric perturbations of BPS monopoles and JZ dyons which give rise to finite angular 
momentum are either unbounded at the origin or at infinity. 



VIII. CONCLUSIONS 



In this paper we have presented a gauge invariant approach to the stationary perturbations of Julia-Zee dyons and 
BPS monopoles. Restricting our attention to axisymmetric perturbations, we have found three sets of modes in each 
parity sector. 

• Electric perturbations: These are manifestly gauge invariant, since the electric background field vanishes after 
a hyperbolic rotation. There exist even parity perturbations with finite angular momentum; however, these are 
either not well-behaved at the origin or at infinity. The same is true of the odd parity perturbations, which give 
rise to axial deformations. 

• Non self-dual magnetic perturbations: These are perturbations which satisfy the linearized field equations, but 
are not at the same time subject to the linearized Bogomol'nyi equations. As the corresponding background 
field vanishes, the non self-dual magnetic perturbations are also gauge invariant. Like in the electric case, there 
exist even parity perturbations with finite angular momentum. However, neither these nor the odd parity modes 
are well behaved. 

• Self-dual magnetic perturbations: These have been investigated before and are known to be physically not 
acceptable, apart from the translational modes. Moreover, general considerations show that self-dual modes 
cannot contribute to the angular momentum. The fact that all solutions to the linearized BPS equations are 
known in closed form is, however, very useful to reconstruct the physical fields 8H and 5 A for the non self-dual 
modes. 
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In conclusion, we would like to emphasize that the distinguished properties of the BPS background are very critical 
to the methods developed in this paper. Whether the main result - the fact that there exist no rotational excitations 
of Julia-Zee dyons and BPS monopoles - generalizes to more general background configurations is an open problem. 
In particular, the effect of a Higgs potential and the coupling to gravity need to be investigated for the excitations of 
Julia-Zee dyons. 

APPENDIX A: BPS BACKGROUND 

The a pprop riate way to treat axial perturbations in gauge theories is by using the isospin harmonics introduced 
in Sect. [VA. It is, therefore, suited to write the background fields in terms of the spherical su(2) basis r r , t&, t v , 



defined by 

T r = t ■ r , dr r = Ttfdi? + t v sint?d</?, (Al) 

where r = r/r is the radial unit direction, and r = <r/(2i). The commutation relations of the Pauli matrices imply 
[r r , Tg] — t v (and cyclic), from which one obtains the formulas 

[r r ,dr r ] = — idr r , [dr r ,*dr r ] = 0, (A2) 

[dr r ,dr r ] = [*dr r ,*dr r ] = 2r r df2, (A3) 

where * denotes the Hodge dual with respect to the standard metric of the two-sphere S 2 . It is helpful to recall that 
the radial unit direction r is a vector valued eigenfunction of the spherical Laplacian with eigenvalue 2, implying 

d*dr r = -2r r dfi. (A4) 

In terms of r r , the "Witten ansatz" for the spherically symmetric connection one- form assumes the simple form 

A = [1 - w(r)] *dr r , (A5) 

since *dr r = r~ 2 (r x t) • dr. Using the commutation relations (jA^), the gauge covariant derivatives of r r , dr r and 
idr r become 

Dr r = w dr r , Ddr r = 0, (A6) 

DSdr r = -2wr r dfl (A7) 



The Bogomol'nyi equations, *F — DH, can easily be written out by using Eqs. (A7), the ansatz H — h(r)r r , and the 
formulas *(dr A *dr r ) = — dr r , *dO = dr/r 2 for the three-dimensional Hodge dual. One finds 

w 2 - 1 

* F = w' dr r H — r r dr, (A8) 

DH = hw dr r + ti r r dr, (A9) 
which yields the well-known first order equations (^5|) for w(r) and h(r). 

APPENDIX B: THE 2 + 1 DECOMPOSITION 

The axial perturbation equations for a static, spherically symmetric su(2) valued function involve the three- 



dimensional gauge covariant Laplacian with respect to the gauge potential (A5). As the latter is tangential to 
S 2 , the three-dimensional gauge covariant derivative operator is 

D = dr A d r + D, (Bl) 

where 
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D = d • +[A, ■ ], with d = dz? A dv + dtp A d v . (B2) 
For an arbitrary Lie algebra valued function / we thus have 

* D * D (i)=l( 92 r + ^* t) ^f, (B3) 

where the factor 1/r is introduced for convenience. (Here we have used *dr — r 2 d£l and *D/ = — dr A *D/.) The 
above formula enables us to immediately write down the 2 + 1 decomposition of the electric perturbation equation 
( ^l| ) . With / = n5$ this becomes 

d r + 4*^*6 ) (r5$) = -[H, [H, r<5$] ] . (B4) 



The 2 + 1 decomposition of the (first order) magnetic equations (22) with respect to the ansatz 

5B = X;bdr + B (B5) 
r j, 

was given in Sect. |VB| ; see Eqs. (fig). (Here B denotes an su(2) valued one-form tangential to S 2 , and b is an su(2) 
valued scalar field.) We owe the proof of the assertion that b is subject to the same second order equation as the 
scalar electric perturbation 5$. In order to see this, one applies *D* on the second, and d r on the third equation in 
(ft8|). A short calculation yields 

dl + -4iDSD ) b = -*[(Di? - *A') , B] - [H, *BB]. 

r z J 

The 2 + 1 decomposition of the Bogomol'nyi equation gives T>H — *A', implying that the first commutator on the 
RHS vanishes. By virtue of the first equation in (Q, the second commutator becomes [H, [H,b]], which yields the 
result 

d 2 r + -i*D*D ] b=-[H, [H, b] } . (B6) 



Hence, the equation (B4) for the scalar electric perturbation, rSQ, coincides with the second order equation ( [Bq ) for 
the scalar part of the magnetic perturbation, b = r 2 (dr,B). 



APPENDIX C: HARMONIC ANALYSIS 



By virtue of the above decompositions, the task of writing out the perturbation equations reduces to the problem 
of computing the gauge covariant derivative D of su(2) valued functions and one- forms over S 2 . We have already 
argued in Sect. IV A that the J = 1 sector is spanned by the three scalar harmonics X, Y and Z, defined in terms 
of r r and K = cosfl [see Eq. (p0|)], and the four one-forms dX = — \A2dy, idX, dZ and *dZ [see Eq. (|32|)]. Instead 
of the latter, it is very convenient to use the linear combinations r r dK, Kdr r and their duals. The entire harmonic 
decomposition is then obtained form the formulas 



DX = T r dK + wKdr r , 
V2BY = -WT r dK-K dr r , 
V2BZ = w T r *dK - Kidr r 

for the covariant derivatives of the scalar basis, and the relations 



(CI) 



it)(T r dK) = -wV2Z, 
*D(KdT r ) = V2Z, 
D(T r idK) = w V2Y - 2X, 
T>{KidT r ) = V2Y -2wX 



(C2) 
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for the covariant derivatives of the basis one- forms. (The equation for DX and Eqs. (C2) are immediate consequences 
of Eq. (A7), while the derivations of the expressions for DY and DZ require slightly more work.) 

As an illustration we compute SDSD(r6$), where we use the expansion @ to write r5<P = cp-X + <j) + Y + 4>Z. 
For the first term we find, for instance 

*D*D (4>-X) = 0_iD [r r *dK + w K*dr r ] 
= 2V2w Y - 2 (w 2 + 1) 

A similar computation for the second and third term gives the result 

JDiD (n5$) = [-2 O 2 + + 2V2wfa 



X 



2\/2w4>- - {w 2 + l)4> + 



(w 2 + 1) 



Y 



which, together with the 2 + 1 decomposition formula (B4) and [H, [H, (rJ$)] 
perturbation equations (Ha) and (Efl). 



(C3) 

-{4>+Y + 4>Z), yields the desired 



APPENDIX D: GAUGE TRANSFORMATIONS 



In this Appendix we show that there exists a gauge for which the perturbations 5H and 8 A assume the expansions 
( p3| ) and (|54j), respectively. We also establish that the coefficients are gauge invariant, up to the residual gauge 
transformations given in Eqs. (^5|) and (^6|). For simplicity, we focus on the even parity sector; the manipulations for 
the odd parity sector are completely analogous. The general expansions for 5H cvcn and 5A avcn are 

6H cvcn = + 7+r, (Dl) 



SA C 



ao Z dr + ot\ T r *dK + a.2 K*dr r , 



(D2) 



where the bars have been introduced to tell the amplitudes apart from the ones introduced in Eqs. 
Under a gauge transformation with an su(2) valued function \ one has 



SH 
8 A 



SH 

5A- 



-[H, X ], 



I and @). 



(D3) 



where, as usual, H is the background Higgs field and D the covariant derivative with respect to the background 
potential A. The strategy is to write 5H and 6A as sums of a pure gauge and an (almost) gauge invariant part. For 
S A this is achieved by a partial integration of the radial part, and by using the expressions ( pl| ) for the covariant 
derivatives of the isospin basis. The radial part of SA even can be written as 



ao Z dr = D 



aodr 



DZ / aodr . 



where we have used the fact that DZ = DZ. (Recall that D — dr A d r + D, and that the isospin harmonics are defined 
over S 2 .) Now using the expression (CI) for DZ brings SA even into the desired form: 



SA C 



-Dx 



w 



1 1 j - —= I a dr 



T r idK 



a < + — / n,,cir 



K*dr r , 



where 

X = -Z J a dr. 

In order to separate a pure gauge term from SH even , we use [r r , Z] = Y and H = hr r to write 



(D4) 



(D5) 
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5H e 



-[H,x]+1-X 



7+ 



agdr 



Y. 



(D6) 



with x according to Eq. (|D5|) . Hence, after a gauge transformation with x, the general perturbations (Dl) and ( |D2| ) 
assume the form (|53|) and (|54|), respectively, where the coefficients are related as follows: 



7- 



7- 



7+ = 7+ 



a dr, 



w f - j 1 / 

ai = ai — — ^= y aodr, q?2 = a.2 + —j= \ n.,ur. 



(D7) 



It is clear form the above reasoning, and not hard to verify, that the amplitudes without bars are gauge invariant, up 
to residual gauge transformations with 



Xo = ciX + c 2 Y + c 3 Z, 



(D8) 



where ci, C2 and C3 are arbitrary constants. Since only the last term is relevant to the even parit y se ctor, we have 
D X c Q vca = c 3 (wT r *dK - K*dr r W2 and [H, 

Xo ven ] = c 3 hY. Using this in the transformation laws (Pq) for the even 
parity perturbations ( |53| ) and (|54|), we conclude that 7_ is gauge invariant, while 7+, oi\ and a 2 transform according to 
Eqs. (pq ) under the residual gauge transformations. A completely analogous reasoning establishes the transformation 
laws ( pq ) for the odd parity sector. 



APPENDIX E: EVEN PARITY ELECTRIC PERTURBATIONS 



In this Appendix we briefly show how the two coupled second order equations ( p5| ) for and 0+ can be trans- 
lated into the inhomogeneous second order equation ( p8| ) for E, defined by S' = h<p'_. The procedure involves two 
integrations. The first integration is achieved by the observation that Eqs. (Ea) can be cast into the form 



V<t>- 



(El) 
(E2) 



where we have introduced the short-hand /1 = w + l/w. Since the RHS of the above equations are equal, and since 
both LHS are exact derivatives [</>"/i — <t>+n" = [fi 2 {4>+ / an integration yields the following first order relation 
between </>_ and <fi+: 



1 

71 



(E3) 



where fc 3 is an integration constant. We now solve Eq. (E2) for perform a derivative and use the result on the 

LHS of Eq. ( E3 ) . This yields the following second order equation for : 



1 



- 2(j)'_ = 4 



ur 



1 



(E4) 



which shows that 0_ = const, is a solution of the system (El), (E2). 

Our aim is to integrate Eq. (E4) once more. In order to see that this is possible, we introduce the variable S 
according to definition (|S7|), and note that the term in front of (<f>'_— k 3 ) can be written in the form [(w 2 — l)/(w 2 + l)]' /h. 
Hence, with 



h(f>'_ 



w 2 - 1 
w 2 + 1 ' 



(E5) 



Eq. (E4) assumes the form 





ft'Y 


hi 


\ h )_ 



(E6) 
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It is not hard to perform the differentiations and to rewrite this third order equation for S in the form 



E + ~ E 



- 2 



2E' - kaa? 



= 0, 



where each of the three pairs is manifestly an exact derivative. Integrating the above expression and multiplying the 
result with h! j a eventually yields 



E"-2^E'-2-E 
h a 



ti 

-2fc kati, 

a 



(E7) 



where fco is a further integration constant. Since a = r 2 h' /(w 2 + 1), this is the desired inhomogeneous second order 
equation ([B8|). We recall that the four parameter family of solutions to Eq. (E7) is 



E = ^ fciE 



(E8) 



where the sum runs from to 3, and where E^ ^ = 1 , E^ 3 ) = — h 2 /2, and E^ 1,2 ) are the two (nontrivial) solutions 
to the homogeneous part of equation (E7). The four independent solutions (|89|)-(|9l|) to the original system jiq ) are 
finally obtained from Eq. (§70. 



APPENDIX F: ELECTRIC CONTRIBUTION TO THE ANGULAR MOMENTUM 



In Sect. IV C we have argued that the total angular momentum can be expressed in terms of the perturbation 
amplitudes at the origin and at infinity. While we have established this result for the magnetic contribution ( |39| ) , we 
still owe the proof of the formula ([l3]) for the electric part ( |38| ) . In order to show that the bracket in the integrand in 
Eq. ( [4l| ) is an exact radial derivative, we first perform a partial integration in both terms, which yields 



rh<p- pic' ' (f 

v2 



- h(r<f>-)' + V2whfa 



(Fl) 



In order to show that the last two terms on the RHS combine to an exact derivative, we use the first perturbation 
equation in (p!5| ) to express cf> + in terms of and <$>"_. Also using the background equations for w and h, we then 
have 

r 2 h 

— h (r<^_) + y2wh(j)+ = — — (f>'!_ — rh <j)'_ + w 2 h cf>— 



■ ^ + [w 2 - 1 + rh] 4>'_ + [2w 2 h + (rh)'] 

( r 2 h \' 
-{—<!>-) +{w 2 -l + rh)4>- 



Using this on the RHS of Eq. (Fll) gives the desired formula 



rh'cf)- 



r i w 



V2 



(1 - w 2 - 2rh) cj>_ + r 2 h(f)'_ + V2wrh(/) + 



(F2) 



which was used in Sect. [V C to establish the result ( 
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